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t:;;^- ' Abstract 

o : 

\ The separation of the Schrodinger equation into a Markovian and an interference 

■ term provides a new insight in the quantum dynamics of classically chaotic systems. 

'"^^1 The competition between these two terms determines the localized or diffusive char- 

^ \ acter of the dynamics. In the case of the Kicked Rotor, we show how the constrained 

^ ■ maximization of the entropy implies exponential localization. 
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. 5t , In the last few decades the field of Quantum Chaos has drawn the attention 

of researchers in several areas of science. Many interesting phenomena were 
studied in systems such as atom traps and microwave cavities. Furthermore, 
the recent advances in technology that allow to construct and almost perfectly 
preserve quantum states, has opened the field of quantum computation, where 
chaotic effects and their control play an essential role (jlj). 

In this work we develop a different and general approach to the subject of 
dynamical localization (DL) and the related issue of quantum diffusion. This 
approach leads to an improved understanding of why DL takes place in some 
systems while in others quantum diffusion continues for ever. The path we 
follow consists in rewriting the Schrodinger equation in a form in which the 
part responsible for DL is separated from the part responsible for quantum 
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diffusion. As we shall see, the terms responsible for quantum diffusion have 
the form of a master equation, typical of Markovian processes, while the other 
part consists of interference terms required to preserve the unitary character 
of the quantum evolution. The master equation has been extensively used by 
W. Zurek and co-workers Q) to study the emergence of a classical behavior 
due to environment-induced decoherence. For the systems considered in the 
present work the environment does not play a role, and the decoherence, if 
present, is self-generated by the dynamics. 

Consider a quantum system described by a generic time- dependent Hamilto- 
nian of the form H{t) = Hq + V{t), where Hq is the time-independent 
part, with known eigenstates and eigenvalues satisfying H()\k) = Ek\k). We 
write the wave function as = J2k (^k(t)\k) . The time- dependent part of 

the Hamiltonian, V{t), induces transitions between the eigenstates of Hq. The 
actual quantum unitary evolution of is given by the Schrodinger equation. 
However, if discrete times tn are considered, then a quantum map is obtained 



afc(t„+i) = J2Uk?(^i{tn)^ (1) 
I 

where f/^"^ = (fc|[/(t„+i, and U{t n+i,tn) is the evolution operator con- 
necting the state at time t„ with the one at time tri+i- In terms of U, the 
occupation probability Pkitn) = \0'k{tn)f can be expressed as 



P,(Wi) = Y: UlTulT(^iitn)a*M. (2) 

l,m 

If we separate the diagonal terms in its rhs, eq. (2) may be written as 



Pfc(t„+l) = ^ TklPlitn) + f3k{tn) , (3) 
I 



with 



A(tn)= E Ul':^Ui2*ai{Qa*M, (4) 



_ r(")rr(")* 

l,m 



where we have defined Ti,i = 



(n) 2 



as the transition probability. This is 
meaningful because T^i > and J2kTki = J2iTki = 1- Thus, if the term Pk 
could be neglected in eq. (3), the time evolution of the occupation probability 
would be described by a Markovian process in which the transition probability 
k ^ l,m a time Atn = tn+i - tn, is given by Tki. 
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We now consider the transition probability per unit time, defined as 

where 6m is the Kronecker deha. A straightforward calculation shows that the 
quantum evolution equation (3), can be written as 



Pk{tn+l) = Pk{tn) + E (^klPl - WlkPk) At„ + A. (6) 

In this form, two qualitatively different terms can be distinguished: one asso- 
ciated with a Markovian process, i.e. a classical-like diffusion, and the other, 
associated to quantum interference effects. This last term preserves the 
unitary character of the evolution. Its contribution depends on the continuous 
or discrete character of the dynamical response spectrum For systems 

with a discrete spectrum, as the periodic kicked rotor, it is of the same order of 
magnitude as the Markovian term, and becomes responsible for the dynamical 
localization found in such systems. A finite time of the order of l/Ao;, where 
Aa; is the average separation in the frequency response spectrum, is required 
in order to resolve the discreteness of the spectrum. At shorter time scales, 
the discreteness of the spectrum has no effect on the dynamics, the Markovian 
approximation holds and the system "mimics" classical chaos. In the case of 
a continuous spectrum, the terms contributing to (3^ in eq. (4) are scattered 
on the complex plane in such a way that their sum becomes negligible com- 
pared to the Markovian term in eq. (6). Then, the unitary evolution is well 
approximated, for arbitrary long times, by a Master equation 

j^=y: (WkiPi - wikPk) (7) 

'^''n l^k 

which results from eq. (6) if Pk is neglected. 

Assuming, for simplicity, that Wk,k-i = Wk,k+i, eq. (7) may be written as a 
diffusion equation, 

dPk ^ 

dtn 2 dk^ ^ ' 

where the diffusion coefficient is 

oo 

D = 2j2Wk,k+il'. (9) 
1=1 

The differential operators in eqs. (7) and (8) represent the discrete derivatives 

dPk/dtrr = [Pk{tn+l)-Pk{tn)]/^U^nAd^Pk/dk^ = [Pk+l{tn) + Pk-l{tn) ' 2Pk{Q] /I 

Note that, according to eq. (8), an initial gaussian distribution evolves as a 
spreading gaussian. 
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As it is well known [6j, a Markovian process described by a Master equation 
such as (7), satisfies the H-theorem, which assures that the thermodynamic 
entropy, S = —J2k Pk^T^ Pk can only increase. On the other hand, in a unitary 
evolution the entropy calculated from the density operator, S= - trace(plnp), 
must be constant. We note that this expression is equivalent to the entropy 5* 
only if the index k can be associated to the eigenstates of the complete Hamil- 
tonian H{t). Obviously in this situation the Master equation does not hold. 
Since in our case we have used the representation of eigenstates of Hq, and not 
of H{t), the thermodynamical entropy 5* corresponds to a "coarse-graining" 
of the entropy S. In the rest of this paper, we refer to the thermodynamic 
entropy S simply as "entropy" . Even though it is possible to define quantum 
dynamical entropies ((71), as long as the Master equation (7) holds it is adequate 
to use the associated thermodynamical entropy S. 

The introduction of the entropy S will allow us to use thermodynamic argu- 
ments to draw interesting conclusions regarding DL. We recall ((8|) that if the 
value of an observable M in state \k) is M^, the maximum value of the entropy 
consistent with a given constraint, 

J2MkPk = {M), (10) 

k 

corresponds to an exponential (canonical) distribution 

= e-^e-^*^^ (11) 

This equilibrium distribution is attained after a diffusive process in which the 
entropy is maximized. The Lagrange multipliers A, Q are determined by the 
constraint (10) and the normalization condition J2kPk = 1, respectively. In 
particular, the condition ^ + (M) = must be satisfied. However, if the 
constraint (10) is removed, the distribution evolves according to a diffusion 
equation for arbitrary long times and the entropy increases without bound, 
provided that the spectrum of Hq is unbounded. 

At this point, it is convenient to consider a concrete example. We choose the 
quantum Kicked Rotor (QKR), one of the first classically chaotic systems 
to be quantum-mechanically investigated (j9|). For this system, the evolution 
operator has the following form (jicl ) 



where k = K/ h is the dimensionless kick strength, Ei = is the kinetic 

energy of the rotor, / its moment of inertia, At„ = tn+i —tn = T the constant 
time interval between kicks and Js is the s^^ order Bessel function. Note that 
in this case, the condition W^^k-i = Wk,k+i is satisfied by (12). Then, as long as 
the term /5fc can be neglected in eq. (6), the evolution is given by the diffusion 
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equation (8). The diffusion coefficient, obtained from eq. (9) using standard 
properties of tlie Bessel functions, is D = k^/2T , wliich is consistent with the 
classical result. 

The average energy of the QKR is given by E(tn) = ^^^=-00 ^kPkitn)- A 
direct calculation based on eq. (6), shows that the energy change across two 
consecutive kicks, AEn = E{tn+i) — E{tn), can be written in terms of the 
transition probabilities WkiAtn as 

00 00 
AEn= E EiPkiQWkiAtn + E 

fc,Z=— 00 Z=— 00 

00 

= E Eif3i{Q. (13) 

l=—oo 

The first term in the rhs of this expression describes the classical diffusive 
increase in the energy while the second term, involving the sum of/3;, accounts 
for quantum interference effects and is calculated explicitly in (p). In Fig. 1, 
we show the contribution of each term in eq. (13) to the energy of the QKR, as 
a function of the number of kicks, n. The dotted line in this figure corresponds 
to the contribution of the interference term, J2i EiPi(tn) and the dashed line 
to the contribution of the Markovian term. This last contribution has a slope 
and is coincident with the classical linear diffusive increase. After ~ 50 
kicks, the cumulative effect of the quantum coherence begins to cancel out the 
diffusive growth of the energy. In eq. (13), At„ = T and the parameters are 
K, = 21.0 and hT/I = 1. The energy is in units of ff/2I and the smoothing 
was obtained by averaging over one hundred initial conditions, corresponding 
to the fifty lowest eigenvalues of Hq. We see that the contribution of the 
interference term is negligible until dynamical localization sets in, so that this 
term may be identified as the one responsible for DL. 

During the quantum diffusion which precedes DL, when the interference term 
is negligible, the quantum evolution mimics a Markovian process, and the en- 
tropy increases, as shown in Fig. 2. However, the discrete quasi-energy spec- 
trum of the QKR Floquet operator (12) introduces a constraint that limits the 
entropy increase. In fact, as is well known ([lol). only a finite number (Lq on the 
average) of eigenstates of Hq are required to describe an eigenstate of the Flo- 
quet operator. If we initialize the system in the eigenstate \ko) of Hq, the sub- 
sequent evolution can only involve eigenstates \k) of Hq such that \k — ko\ < L. 
This fact represents a constraint on the entropy increase, of the type given by 
eq. (10). The actual value of L depends on the initial condition, but it must 
obey (L) = Lq, where the average is over different initial states |/co)- For each 
ko, we take L, i.e. the maximum distance \k — ko\ in angular momentum space 
present in the dynamical evolution, as the observable M of eq. (10). Then, the 
corresponding constraint is (L) = {\k — ^olmax) ~ -^o- From this constraint 
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Fig. 1. Contribution of eacii of the terms of eq. (13) to the average energy of the 
kicked rotor. See text for details. 

and the normalization condition for the occupation probability, we obtain the 
Lagrange multipliers A = sinh~^(l/Lo), f2 = In Lq + (Lq + 1)^/^ . Then, the 
probability distribution (11) resulting from the constrained maximization of 
the entropy is 



1 Ifc-fcpl 



(14) 



Lo + [Ll + 1)1/2 2Lo 

because Lq ^ 1. This exponential profile for the localized wave function has 
been verified both numerically ^1^ and experimentally (11) for the case of the 
QKR. When the maximum value of the entropy consistent with this constrain 
is attained, the interference term in the rhs of eq. (6) becomes non-negligible 
and so the master equation approximation does not hold. Then the entropy 
remains constant (see Fig. 2). 

The origin of DL may thus be traced to the discrete nature of the dynamical 
response frequency spectrum (or equivalently for the QKR, the quasi-energy 
spectrum) , which determines the degree of randomness of the values in the 
complex plane. This is discussed in further detail in (j^. 

To clarify the relation of DL to the behavior of the second term in the rhs of 
eq. (13) we now consider a simple modification of the QKR, in which the time 
intervals At„, are randomly chosen. In this case eqs. (12) and (13) are valid, 
but the contribution of the interference term in the rhs of eq. (6) is negligible 
for arbitrarily long times since the dynamical response frequency spectrum is 
continuous and this system does not show DL (|5j). In this case, the evolution 
is always Markovian and it can be described by the diffusion equation (8) with 
the diffusion coefficient given by eq. (9). 
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Fig. 2. Entropy S = — P^lnP^ as a function of the kick number (log scale) for 
the quantum kicked rotor with a fixed time interval between kicks (thick line) and a 
random time interval between kicks (thin line). The parameters and averaging over 
initial conditions are the same as in Fig. 1. 

One might think that these diffusive dynamics is an artifact resulting from 
the randomness of the time intervals between successive kicks. That this is 
not the case may be seen by considering the double Kicked Rotor (QDKR), 
a quasiperiodic version of the QKR. The QDKR is obtained when a second 
series of pulses is applied to the rotor, with an irrational ratio between the 
periods of both trains of pulses. The evolution operator (12) also holds for the 
QDKR, but the time intervals At„ now form a pseudo-random sequence. The 
response spectrum of this system is dense and it undergoes unlimited quantum 
diffusion at a rate consistent with the classical diffusion rate i). This implies 
that the contribution of the interference term in eq. (6) is negligible also in this 
case. In both cases, the probability distribution Pkitn) evolves as a spreading 
gaussian. The long-term evolution of the entropy is markedly different for 
these systems than for the periodic case. In Fig. 2, the time evolution of the 
entropy S for the QKR and the random kicked rotor are compared. In the 
case of the QKR, the entropy stops growing when the constraint is satisfied, 
while it grows without bound in the other case. 

Finally, we remark that the approach of separating the Schrodinger equa- 
tion into a master equation part supplemented by a term which takes into 
account quantum interference effects, provides a new insight on quantum dif- 
fusion and dynamical localization. As we have shown, the localized or diffusive 
character of the dynamics can be understood as the result of the competition 
between these terms of the evolution equation. The Markovian term domi- 
nates the dynamics for a time which depends on the topology of the dynam- 
ical response spectrum. The general procedure has been illustrated using the 
quantum kicked rotor, in its periodic, quasi-periodic and random kick ver- 



7 



sions. In particular, we have explained the exponentially localized probability 
distribution of the QKR using simple, well accepted results from the theory 
of stochastic processes. This exponential distribution results from the con- 
strained maximization of the coarse-grained entropy. The localization-length 
constraint is due to the discrete nature of the dynamical response spectrum of 
the QKR. In the other cases considered the dynamical response spectrum is 
dense, the localization length is infinite {i.e. the constraint is absent) and the 
entropy increase continues indefinitely, as is characteristic of a diffusive pro- 
cess. We have shown, that in this case the Schrodinger equation is equivalent 
to a diffusion equation with the classical diffusion coefficient. 

The approach that wc have presented here may be used to describe the dy- 
namics of an arbitrary quantum system in which there is quantum diffusion. 
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